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Relaxed quaternionic Gabor expansions at critical 

density 

Stefan Hartmann 


Abstract 

Shifted and modulated Gaussian functions play a vital role in the representation of signals. We extend 
the theory into a quaternionic setting, using two exponential kernels with two complex numbers. As a 
final result, we show that every continuous and quaternion-valued signal / in the Wiener space can be 
expanded into a unique series on a lattice at critical density 1, provided one more point is added in 
the middle of a cell. We call that a relaxed Gabor expansion. 


1 Introduction 

Phase, as have shown A. V. Oppenheim and J. S. Lim in the 1980s, codes all relevant details of an image 
[16]. However, reasonable phase terms for higher dimensional settings, e.g. for 3-D image stacks or for 
6 -I- 1 D (3D color films, e.g. from confocal microscopy) have long been an open question. Among the 
different concepts in 1999 T. Biilow introduced the hypercomplex signal in [6]. The concept of hypercomplex 
signal was afterwards widely used in image processing and even applied to movies (see for instance the 
paper of W. Chan, H. Choi and R. Baraniuk [8, 9]). Billow’s approach is based on a specific type of a 
quaternionic Fourier transform. There exist several possibilities how one can define a quaternionic Fourier 
transform, see e.g. T. A. Ells article [10]. Among these, the right sided version has been already widely 
discussed by E. Hitzer, M. Bahri et al. ]1, 2, 3] and others. Billow’s concept is based on the two sided 
version, where the signal is wedged between two exponential kernels. This idea leads to more symmetric 
properties [11]. For the mathematical background we refer to [5, 7]. But for practical application using 
only a Fourier expansion there is a major drawback. The quaternionic Fourier transform lacks an important 
property, this is the locality of phase behavior. In order to overcome this issue we multiply the signal by a 
real-valued window function with compact support and then apply the quaternionic Fourier transform. That 
way we obtain a representation both in space and phase and the so-called windowed quaternionic Fourier 
transform (WQFT). But in applications we can only work with a discrete set in the parameter space such 
as a lattice. A discretization of the WQFT leads into Gabor frame theory. An important question in that 
theory is, how dense we have to sense the signal, or in other words, which is the minimum size a cell in the 
lattice should have. In the non-quaternionic theory it is well-known, that under certain requirements to the 
window function, the area of a cell has to be smaller or equal than 1 for a perfect reconstruction, whereas 
the case of greater than 1 never leads into a perfect reconstruction. For the special case of the Gaussian 
as window function, the case equal to 1, also referred to as critical density is of main interest, since the set 
of coefficients is complete, but does not constitute a Gabor frame. V. Palamodov showed that there exist 
a unique expansion into shifted and modulated Gaussian windows for the critical density, provided one 
point is added in the middle of a cell [17]. Unfortunately, the same can not be said about the quaternionic 
case. Up to now there is no result regarding the critical density. The principal goal of this work is to 
construct Gabor expansions with critical density, so-called relaxed Gabor expansions, which is essential for 
the practical application of quaternionic Gabor expansions in Image Processing. 

The paper is organized as follows: In Section 2 we briefly review some results on quaternionic algebra, the 
used function spaces and the WQFT. In Section 3 we introduce Gabor methods, such as the Gabor transform 
and the Gabor series. The fourth section treats the quaternionic Zak transform, a natural extension of the 
Zak transform, which is the main tool for the results. In the last section we show that a quaternionic signal 
has a unique P' expansion at critical density, when one more point is added to a lattice. Furthermore, a 
formula for the coefficients of the Gabor expansion is obtained. 
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2 Preliminaries 

2.1 Quaternions 

Quaternions are a non-commutative extension of the complex numbers into four dimensions. The quaternion 
algebra HI is given by 


M = {q\q = qo + iqi+ jq 2 + kqs, with qo, gi, 92 ,93 G R}, 
where the elements i,j,k satisfy 


ij = —ji = k, = ijk = — 1. 

We can write a quaternion q also as a sum of a scalar Sc[( 7 ] = 90 G and a three-dimensional vector 
Vec[ 9 ] = iqi + jq 2 + kq^, which is often called pure quaternion. The conjugate is an automorphism on H, 
given by 9 I—>■ 9 = Sc[ 9 ] — Vec[ 9 ]. The modulus I 9 I of a quaternion is defined by 

kl = = \j Qo + <li + <l2 + <?3- (1) 

It also holds Sc| 9 | < kj and 99 = \q\^. The inverse of a non-zero quaternion 9 is given by The non¬ 
commutativity of the quaternions means that in general p~^q ^ such that ^ in general does not make 

sense. Nevertheless, for an easier reading we will use the latter notation from time to time but we will make 
clear what is meant by it in each case. The following properties also should be brought up 

p^=q-p and |p ■ 9 ^ = Ipp • kp. 

Although quaternions are non commutative, the cyclic multiplication has to be mentioned, since it is going 
to be a useful tool. 

Lemma 1 (Cyclic Multiplication [14]) 

For all 9 , r, s G H holds 

Soke’s] = Sc[rs9] = Sc[s9r]. (2) 


2.2 Function spaces 

In the following we use the space L^(]R^,]HI), an immediate generalization of the Hilbert space of all square- 
integrable functions, which consists of all quaternion-valued functions / : —>■ HI with finite norm 

II/II2 = (^^J/(x)kd^x^ <00 

where d^x = dxi da ;2 represents the usual Lebesgue measure in K.^. The L^-norm is induced by the 
symmetric real scalar product 

{f,9) = lf ( 5 (x)/(x)-k/(x)g(x))d^x = Sc / /(x)g(x)d^x (3) 

which makes the space a real linear space. 

Remark 1 ([15]) 

It is also possible to define a quaternion-valued inner product on L^(]R^,]HI) with 

if, 9 )= f /(x)g(x)d2x 

JR2 

in which case one obtains a left Hilbert module. One has to be careful with the properties, since for p G H 
holds 

ipf, 9 ) = Pif, 9 ) but (/, gp) = if, g)p. 

Both inner products lead to the same norm. 
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We denote the two-dimensional cube [0,1]^ briefly by Q. Moreover we need the notion of a Wiener-Amalgam 
space in the quaternionic setting. The Wiener-Amalgam space comes from the periodization of functions, 
the main idea for the definition arises from a joint categorization of local and global behavior. 

Definition 1 (Wiener-Amalgam space [13]) 

A function / G L°°(]R^,]HI) belongs to the Wiener-Amalgam Space W = W(K.^,L°°) if 

ll/lk= E Wf-TnXQW^ 

neZ'i 


is finite. 

The subspace of continuous functions is denoted by Wq. 


2.3 The windowed quaternionic Fourier transfom 

For the analysis of signals, especially for results on the phase of an image, the windowed quaternionic Fourier 
transform will be our main tool. It is useful to utilize two exponential factors, which can be put either on 
the same side of the signal (one sided) or they can sandwich the signal (two sided). Since we are interested 
in applications using the hypercomplex signal, we use the two sided version which has some more symmetric 
features than the one sided one. The reader should note that these two Fourier transforms differ, owing to 
the non-commutativity, when a quaternion-valued signal is analyzed. 

Definition 2 (Two sided quaternionic Fourier transform (QFT) [11]) 

The two sided quaternionic Fourier transform of / G L^(R^,IHI) is the function iFqif) : ^ HI defined by 

iFq{f){uj) = f{uj) = / exp(-27ria;iWi)/(x)exp(-27rja:2W2)d^x. 

Remark 2 ([11]) 

We can reconstruct the signal / from its QFT with 


/(x) = J" ^(J,/)(x) = / exp(27ria;iWi)/(u;)exp(27rja:2W2)d^w. (4) 

There exists also a Parseval Theorem for the QFT, which is given in the following Theorem. 

Theorem 1 (Parseval Theorem [11]) 

For /, (/ G (R.^, H) holds 

{f,g) = 

Proof. We use the reconstruction formula (4) and the cyclic multiplication property (2) 


(/, g) = Sc [ 

/(x) 5 (x) d^x 



= Sc / 

/ exp(27rza:ia;i)/(a;) exp(27rjx2W2) g{x) d^x 


dR2 

= Sc / 

/(w) / exp(27rja;2W2)(/(x) exp(27ria:iWi) d^w d^x 

dR2 

JR2 

= Sc / 

/(w) / exp(—27rza;ia;i)(7(x) exp(—27rja;2W2) d^x d' 

dR2 

dR2 

= Sc / 

f{uj)g{uj) d^uj 

dR2 


= (/.?)■ 



The cyclic multiplication is essential for the proof. So far there is no Parseval Theorem for the quaternion¬ 
valued inner product. As a simple corollary we obtain the Plancherel formula, which states that the energy 
of a signal is preserved under the quaternionic Fourier transform 

Il/ll2 = ll7lli (5) 
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Nevertheless it is in general not possible to obtain local information of the frequency behavior of the signal 
/ only by looking at /. In order to obtain local information of the frequency behavior, there are several 
possibilities. One could be to cut the signal into pieces and then apply the QFT. Another idea is to multiply 
the signal with a real-valued function of compact support and then apply the QFT. The hrst method leads 
into the Gibbs-phenomena, therefore we use the latter. 

Definition 3 (Windowed quaternionic Fourier transform [11]) 

The windowed quaternionic Fourier transform (WQFT) of / € L^(K.^,]HI) with respect to a non-zero real¬ 
valued window function g G L^(]R^,IR) is defined as 


Qgf{h,uj)= / exp(—27rta;iu;i)/(x)5(x — b) exp(—27rja;2W2) d^x. 


/R2 


( 6 ) 


We can reconstruct the signal from its WQFT. 

Proposition 1 (Reconstruction formula [11]) 

Let g € L^(K.^,K.) be a non-zero real valued window function. Then every / G L^(R^,]HI) can be fully 
reconstructed by 


/(x) 




exp{2TTixiOJi)Qgf{h, oj)g{x - 


b) exp(27rja:2a;2) d^w d^b. 


(7) 


By looking at the definition of the WQFT (6), we can see that the signal / stays untouched. The translation 
and modulation work on the window g, which can be interpreted as a T^-kernel. This gives the first 
motivation to rewrite the WQFT with the help of an inner product. The problem is that we have to find a 
way to express the kernel. In order to do that, we use the concept of carriers, introduced by Shapiro et al. 
[ 18 ]. 

Definition 4 (Carrier [12, 18]) 

We define for two quaternions p,q gM. the right Cr and left Ci carrier operator 


Cr{p)q = qp and qCi{p) = pq. 

Lemma 2 (Properties of the carrier) 

We have for the carriers the following properties with p G H 
(a) Cr{p) = Clip) and Clip) = Crip), 
ih) Crip)l = ICiip) = p. 


Now that we are equipped with the necessary tools we can generalize the translation and modulation for our 
setting. 

Definition 5 

For b, w G we define the following operators 

Tbg(x) = 5(x - b) 


and 


Mc^gix) = exp(27rjw2a^2)ff(x)C'r(exp(27rfa;ia;i)). 


This allows us to write the WQFT in the following manner 


Qgfih,co) = if,M^Tkg). 

For the rest of this paper we employ the two-dimensional Gaussian function as window function, since it 
minimizes the uncertainty principle, i.e. it has the best localization in space and frequency. For shorter 
notation we use 


eA(x) = M^Tb exp(—7r|xp) = exp(27rja;2W2) exp(—7r|x — bp)Cr(exp(27rza:iWi)) 


and 

e)^(x) = exp(—7r|xp) = C'i(exp(27rza;ia;i)) exp(—7r|x — bp) exp(27rja;2W2) 

with A = (b,a;) being a four-dimensional point in x R^. 
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3 Gabor methods 


In this section we will discuss necessary properties of the Gabor transform introduced in the previous section. 
Let us start with the continuous Gabor transform. 


3.1 Gabor transform 

With the reconstruction formula (7) for the WQFT, we get for the special case g(x) = exp(—7r|xp) the 
following expression 

/(x) = [ [ exp(27ria:iWi)(/,M^rb5)g(x-b)exp(27rja;2W2)d^u;d^b= [ [ (/,eA)e^dA. 

Jr2 7r2 7r2 7r2 


Since every quaternionic function can be split into a sum of real functions, we get the explicit formula: 


/(x) 



exp(27ria:ia;i)(/, e^) exp(—7r|x 


bp) exp(27rja:2W2) d^w d^b 


[ exp(27ria;iwi) [(/, ex) + ex) + {-jf, ex)j + i{-kf, ex)j] 


R2 

exp(—7r|x — bp) exp(27rja;2W2) d^w d^b. 


This formula helps us in the next lemma. For / G L^(IR^,IHI) the function A i—>■ (/,©>) defined in x is 

called Gabor transform. With the next Lemma we show that the Gabor transform is unitary. 

Lemma 3 

For an arbitrary function / G L^(R^,IHI) holds 


[ [ |(/,e,)pd^bd2a; = i||/|p 

JR2 Jr2 4 


and 


/ / / exp(—27ria;iWi)/(x) exp(—7r|x — bp) exp(—27rja:2W2) d^ 

1r2 7r2 JR2 

Proof. Since 


d^bd^o 


|(/,e,)pd2A=||/|p. 


f \{f,ex)\^d^u; 

1r2 

= / Sc / exp(27ria;ia;i)/(x) exp(—7r|x — bp) exp(27rja;2W2) d^x 

JR2 JR2 

= [ Sc|J-,(/b)(-w)P d^w 

= [ Sc|/b(w)P d^x 

JR2 


d^oj 


( 8 ) 

(9) 


with /b(x) = /(x) exp(— 7r|x — bp). We used Plancherel Theorem (5) in (8). Integrating (9) over b and 
changing variables y = x — b leads to 


f Sc f |/(x) exp(—7r|x — bp)p d^xd^b 

JR2 JR2 

= [ Sc / |/(x)exp(-7r|yp)pd^xdV 

JR2 JR2 

= [ |exp(-7r|yp)pdV Sc [ |/(x)pd^x 

JR2 JR2 

= \\\f\\l 
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Finally we obtain 

if, e\) = if, ex) + i{-if, ex) + {-jf, ex)j + i{-kf, ex)j 

to find 

[ [ if,ex)d^hd^co=\\f\\l. 

JB2 JR2 

Now we have to take a look at the discrete version of the Gabor transform. 


□ 


3.2 Gabor series 

In applications we have to replace the integrals in (7) by sums, which results in a discretization that leads to 
lattices in the space-frequency-parameter. In the following we define a lattice where each cell has the volume 
1. We take some elements a,/3 G M’*' such that a/3 = 1 and consider the lattice A C x generated by a 
and (3 

A = {A = am + ,5n, m, n G 7?}. 

Introducing a coordinate b G such that b(a) = 1, then we have for the dual coordinate w(/3) = 1, w G R^ 
in the frequency plane. This definition gives a lattice where the volume of a cell equals 1, the so-called 
critical density. The next proposition shows, that it is possible to reconstruct a quaternionic-valued signal 
with the help of shifted and modulated Gaussian functions in that lattice. 

Proposition 2 

For an arbitrary quaternion-valued sequence {ca} G ^^(A) with A = (b,w) the series 


/(x) = ^ exp(27ria;iWi)cA exp(-7r|x - bp) exp(27rja;2W2) = ^ CAe';^(x) 
AgA AgA 


( 10 ) 


converges in L 2 and the following inequality holds 


1/2 


1/2 


I2 < 0-0 


H |CA, 

AAj=(b2,W2) 


H |CA. 


where exp (-^) • 

Proof. For A = (b, oj) and /x = (d, ry) we have 


A,pGA 


( 11 ) 


= Sc / CAeA(x) c^e'^ix) d^x 
aga mga 

= Sc / 2, exp(27ria;iu;i)cA exp(—7r|x — bp) exp(27rja;2W2) 

aga 

'yP exp(—27rja;2?72) exp(—7r|x — dp)c)7exp(—27ria;ir7i) d^x 
mca 

= Sc / yP exp(27rta;ia;i)cA exp(—7r|x — bp) exp(27rja;2(w2 — 772)) exp(—7r|x — dp)c)7exp(—27ria;i?7i) d^x. 

(12) 
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Computing the a;2-integral of ( 12 ) gives 

/ exp(— 7 r|a :2 — 62 !^) exp(27rja;2W2) exp(—27rja:2??2) exp(— 7 r|a ;2 — ^ 2 ^) dx 2 
Jr 

= / exp(—27ra;2 + 2 Trb 2 X 2 — Tr&f + 2 Trd 2 X 2 — + 2 TTjx 2 {oJ 2 — di)) dx 2 

Jr 

= exp(-7r(&| + ^ 2 )) / exp(—27r(x| — a; 2(&2 + ^2 + W 2 j — r72j))) da ;2 

Jr 

( ri .2 , J2^ , + d 2 + U} 2 j - mjf 

= exp(- 7 r (&2 + d^) H- - -) 

= exp(- 7 rj (62 + d 2 ){uj 2 - m) - ^((^2 - ^ 2 )^ + (w 2 - mY)) 

= exp(- 7 rj (62 + d 2 ){uj 2 - m) - \^j - MjP) 

= exp(- 7 rj(fc 2 - d 2 ){uj 2 - m) - |Aj - d-j?)- (13) 

The last line ( 13 ) holds true, since we have chosen our lattice in a manner that we can use periodicity 
conditions. We insert this in ( 12 ) and get with the help of the cyclic multiplication rule ( 2 ): 

Sc / ^ exp( 27 ria;iu;i)cAexp(- 7 r(a;i - 61)^) exp(-7rjX&2 - d2)(w2 - ?72) - ^((&2 - ^2)^ + (w2 - ?72)^)) 

exp(—7r(a;i — di)^)c^exp{— 2 Trixirii) dxi 
= Sc / ^ CAexp(-7rj'(&2 - d2)(w2 - ??2) - ^((^2 - ^2)^ + (w2 - ??2)^))^ 

exp(—7r(a;i — 61)^) exp(— 7 r( 27 ria:i(u;i — ?7i))) exp(—7r(xi — c?i)^) dxi 
= Sc ^ CAexp(-7rjX&2 - d2)(w2 - m) - Ziih - ^2)^ + {0J2 - ^72)^))^ 


exp(-7rt(6i + di){uji - di) - -^({h - dif + (w2 - d2f))- 

For the estimate we incorporate the discrete Young inequality for a convolution with r = p = 2 and q = 1: 

Y ^nbm-n ^ f XI X 

and the functions 

gx^ = exp(7rj52W2 - ^(&2 + ^ 2 )) 9\i = exp(7ri6iu;i - ^(&? + uif)) (15) 

for Xj = (62,^2) G K X M and \i = (61, wi) G M x K. Moreover, we have 

X = X I^yI = = ^ 0 - ( 16 ) 
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From this we obtain for the norm estimate of / 


- I X! exp(-7rj(62 - d2){uj2 - m) - - ^ 2 )^ + (w2 - V2)‘^))c^j 


x,tJ.eA 
C/ii 


exp(-7rt(&i - fii)(wi - 771 ) - -((61 - di)^ + (w 2 - 


^ XI |cAj exp(- 7 rj (62 - fi 2 )(a ;2 - 772 ) --((62 - ^ 2 )^ + (w 2 - ? 72 )^))c 

A,/^GA 


I 


jcA^c^i exp(-7ri(5i - di){uii - 771 ) - -((5i - di)^ + (wi - 77i)^))| 

XEE cAj 9 \j flEE CA*C^i5Ai-Mi|^ 

f^3 

<Ek.i^E 



2 

2 





Ai Ai 



= EkpiEki Ei^.7pEi^A.nEi^^Aj El 


|CA 

Ai \ A 


Ei^a^P e 


A^j Ai 

|CA.P^ 


(17) 

(18) 


In (17) we use the convolution estimate (14). For the estimate of (18) we use equation (16). 
All this leads into 

//2 / ^ 1/2 

|2 


^-oME 


CA, 


E 


CAi 


□ 


4 Quaternionic Zak transform 

The quaternionic Zak transform is the natural extension of the Zak transform into the two-sided quaternionic 
setting. The classic Zak transform is useful in the analysis of the Gabor transform, of course, the extension 
keeps similar properties [12, 13]. In harmonic analysis the Zak transform is also known as Weil-Brezin map. 
For our results we will invoke the quaternionic Zak transform several times. 

Definition 6 (Quaternionic Zak transform) 

The quaternionic Zak transform of a quaternion-valued function / G T^(IR.^,]HI) is defined as 

^<?/(x,w) = E exp(27riTOia;i)/(x — m) exp(27rj7T7.2W2) x, w G 


In what follows we show some interesting properties of the quaternionic Zak transform. 

Proposition 3 (Quasiperiodicity conditions) 

The transform Zq fulfills 

Zgf{x,uj + 1) = Zq/(x,a;), Zqf{yi+ l,w) = exp(27riwi)Zq/(x, w) exp(27rju;2). 

The next remark concerns the zeros of the quaternionic Zak transform of a modulated and translated 
Gaussian window. 



Remark 3 

We have for ^ = (p, 0)-. 

Z,e^(x,w) 

= exp( 27 riTOia;i) exp( 27 ri 0 i(xi — mi)) exp(— 7 r(a;i — mi — pi)^) 

mGZ2 

exp(— 7 r(a ;2 — m 2 — ^ 2 )^) exp(27rj02(a;2 — m2))(exp(27rjm2W2)) 

= exp(27ri0ia:i) ^ exp(27ri(mia;i — 0 imi)) exp(— 7 r(a ;2 — m 2 — ^ 2 )^) 

exp(— 7 r(xi — mi — pi)^) exp( 27 rj(w 2 W 2 — 02m2))j exp( 27 rj 6 < 2 a; 2 ) 

= exp( 27 ri 0 ia;i) ^ exp( 27 ri(miWi + piWi — ^im'i — Oipi)) exp(— 7 r(a;i + 

m'eZ2 

exp(— 7 r(a ;2 + ^ 2 )^) exp(27rj(TO2W2 + P 2 W 2 — 92m'2 — 6 * 2 ^ 2 ))^ exp(27rj02a:2) 

= exp( 27 rz( 0 ia;i + piWi)) exp(—Trxi) ^ exp( 27 rzmia;i — 27 ra;imi — zrm^^) 

m'GZ^ 

exp(27rjm2a;2 — 2TTX2m2 — exp(— 7 ra; 2 ) exp(27rj(022^2 +P 2 W 2 )) 

= exp(27rz(0ia;i + piWi)) exp(-7ra:i)0j(wi +za;i) 0 j(u ;2 + ja; 2 ) exp(— ttcc^) exp(27rj(02a;2 +P 2 W 2 )) 
with the two Theta series 

0 i(' 2 i) = exp( 27 rimz — zrm^) and Qj{z2) = exp( 27 rjmz — zrm^) with zi G Ci,Z2 G Cj. 

m^Z mGZ 

For the involved Theta series we have the following properties. 

Remark 4 

The Theta series are holomorphic, have a simple zero in the cube Q only for zi = and Z2 = 1 / 2 +j 72 , 

respectively [4]. Moreover, the Theta series satisfy 

0721 + 1 ) = 0 *( 2 i), 0 *( 2 i +i) = exp( 7 r - 27 rzzi) 07 zi) 


and 

0j(z2 + 1) = Qj{z 2 ), '^j{z 2 +j) = exp(7r - 2 njz2)Qj{z2). 

We also have the following continuity property of the quaternionic Zak transform. 

Lemma 4 

For any / G VFo the quaternionic Zak transform is continuous. 

Proof. Given e > 0, there exists a iV > 0 such that 

\k\>N 

Then, the main term exp(27rzfciu;i)/(x — k) exp(27rjA:2W2) is uniformly continuous on compact sets 

of and there exists a <5 > 0 such that 


E 

|k|<Ar 
whenever |x — 


exp(27riA:ia;i)/(x — k) exp{2Trjk20J2) — E exp(27rzfci^i)/(y - k) exp(27rjA:26) 



|k|<W I 

y| + |w — ^1 < S. As a consequence, \Zqf{x,uj) — Z,/(y, ^)| < e and Zqf is continuous. 


□ 
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Furthermore, the quaternionic Zak transform is a unitary mapping. 

Lemma 5 

The quaternionic Zak transform is a unitary operator mapping L^(]R^,EI) onto L’^{Q x QjH). 

Proof. We consider the functions Fin(x, w) = exp(27rimiu;i)/(x — m)(exp(27rjTO2W2)) with m G 1?. These 
functions belong to x QjH)- 

ll'^in||z,2(QxQ) = (-fmjTm) 


= Sc / / exp(27riTOia;i)/(x — m) exp(27rjm2W2) 
mez2 “'<3 “''3 

exp(—27rjTO2W2)/(x — m) exp(—27riTOia;i) d^x d^w 
= Sc / / exp(27riTOiWi)/(x — m) exp(27rjm2W2) 

mGZ2 

exp(—27rjTO2W2)/(x — m) exp(—27riTOia;i) d^x d^w 
|/(x-m)pd^xd^w 


= E 

mGZ2 “'Q “'Q 


For m 7 ^ n, we have with the cyclic multiplication (2) 


{Fm,Fn) 

m,nGZ2 


^ Sc / / exp(27rimiu;i)/(x — m) exp(27rjTO2W2) 


exp(—27rjn2W2)/(x — n) exp(—27riniu;i) d^xd^w 
Sc / / exp(27rimia;i)/(x — m) exp(27rj(TO2 — n 2 )a; 2 )/(x — n) 

exp(—27rzniWi) d^xd^w 

Sc / / /(x — m) exp(27rj(m2 — n 2 )a; 2 )/(x — n) exp(—27rmia;i) 

Jo Jo 


/Q JQ 

exp(27rzmia;i) d^xd^w 


By combining the above results we get that ZmGZ^ Fm converges in L^{Q x Q,!!!) and 

II ^ ^m||i2(QxQ) = X/ ll'^mlli=^(QxQ) = II/II 2 

Thus Zq is an isometry from T^(R^,]HI) to L‘^{Q x (5,]HI). For the rest of this proof we use the Gabor 
orthonormal basis 

TVn.n = exp(27rzmia;i)T’(x — n) exp(27rjm2a;2) 

where X denotes the characteristic function of Q, defined by X{x) = 1 for x G Q and X{x.) = 0 otherwise. 
Direct calculation tells us that 

'2gT’m,n(x, w) = exp(27rz(TOia;i — nia;i))ZqT’exp(27rj(TO2a^2 — ^ 2 ^ 02 )) 

and ZqX = 1 with (x, w) G Q x Q. Since -ZqT’m.n forms an orthonormal basis with the scalar inner product 
(3), Zq maps an orthonormal basis of to an orthonormal basis of L^{Q x Q,]HI). Thus Zq is 

unitary. □ 

The next lemma shows that it is possible to reconstruct the signal / only by the values of Zqf on the cube Q. 
This means we only have to analyze the Zak transform of the signal in the cube to reach global information 
about /. 
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Lemma 6 

If / e 


/(x) = [ Zgf(x,u})d^u} 

JQ 

Proof. We use the definition of the Zak transform and obtain 
[ 2gf{x,u))d^uj 

Jq 

= / exp(27ria;iTOi)/(x — m) exp(27rjw2m2) d^w 

meZ2 

= / /(x) d^w + / 6xp(27riwimi)/(x — m) exp(27rja;2TO2)) d^w 
JQ JQ _in 


m^O 


The first integral is just /(x). In order to calculate the second integral we interchange the order of integration 
and summation and get 

/ exp(27ria;imi)/(x — m) exp(27rja;2TO2)) d^w 

m #0 

= >, / exp(27ria;iTOi) dwi /(x — m) / exp(27rja;2TO2) da ;2 
m^0-^[0d] AoA] 


= 0 . 

This completes our proof. 


□ 


5 Relaxed Expansion 

We show now that an arbitrary function / € Wq can be expanded in a unique Gabor series as in (10) when 
one element in the middle of a cell is added. Denote by # = (1/2,1/2) a point in the middle of the cell. We 
call the set A# = A U {#} the relaxed lattice. We use similar concepts as in [17]. 

First of all, we need to define the sharp Poisson functional, which will be a necessary tool for the results. 

Definition 7 (Sharp Poisson functional) 

The sharp Poisson functional is the series 

l*{f) = (*0.(O))-iZ,/(#)(j0,(O))-i = (19) 

The sharp Poisson functional has the following important property, which is needed to find the coefficient 
at the sharp point ff. 

Lemma 7 

We have 7 ^(e'^) = 1 and 7 ^(e)^) = 0 for arbitrary A G A. 

Proof. The use of the Euler formula for quaternions gives for 7 '^(e(^) 

(*0*(O)) l/2)(0j(O)j) 

= (f0i(O))“^ exp(7rfTOi) exp(—Trm^) exp(— + mmi) exp(7rjm2) exp(— ttto^) exp(^ + 7rjm2)(0i(O)j)“^ 

= (.e.(o))-‘ y; exp(—7rm2)(cos(7r/2) + i sin(7r/2))(cos(7r/2) + j sin(7r/2)) exp(—7rm2)(0j(O)j) ^ 

= (*0.(O))-i(*0.(O))(0,(O)j)(0,(O)jr' 

= 1 
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And for A = (b, w) 

= (* 0 .(O))-iZ,e';,(l/ 2 , l/2)(0,(O)jr' 

= (i0i(O))“^ exp(7ri(&i + wi)) exp(-7r/4)0i(l/2 + i/2)0j(l/2 + j 72) exp( 7 rj (62 + W2))O'0j(O))“^ 

= 0 


since 0^ and Qj are zero at the sharp point, see Remark 4. □ 

In the following theorem we show that the Sharp Poisson functional is well-defined. 

Theorem 2 

For any / € IFo there exists a family of continuous functionals A G A in the space IFo, such that 

./W = ^ exp(27ria;iWi)7#(/)exp(-7r|x-b7)exp(27rja;2W2) = ^ (20) 

xeA* xeA* 

Proof. By Lemma 4, is well-defined in Wo(R^)- Set 

/#(x) = /(x) - exp(-7r(xi - exp(7ria;iWi)y^(/) exp(-7r(x2 - exp(7rja;2W2) (21) 

and define F(x, w) = Zqf^/ZqG^ = exp(7ra;i)(0i(wi-l-ia:i))“^Zq/^(x, w)(0j(a;2-|-ja:2)) exp(7rx2). By Lemma 
8, F is square integrable and has the following double Fourier series expansion 

F(x,w) = ^ exp(27r7a;ipi + u;i?7i))c^exp(27rj'(a;2P2 + W2?72)) (22) 

M=(p.’))6A 

By Proposition 2, the series ^ exp(27ria;iwi)c^ exp(—7r|x — bp) exp(27rja:2W2) = X) converges to a 
function g. 

On the other hand we have with Remark 3 


Zqf^(x,Uj) 

= exp(—7rxi)0i(a;i + ixi) 

M=(’).p) 6 A 


exp(27ri(xipi + 0Jir]i))c^ exp{2TTj{x2P2 + W2??2)) 

M=('/.p)eA 


<dj{uj2 + jx2)exp{—TTxl) 


Which means by the unitarity of Zg that = g, set y^(/) = for /x G A. □ 

Lemma 8 

We have F(x, w) G L'^{Q x Q). 


Proof. The function is an element of IFo- By Lemma 4, Zqf^{x,u)) is continuous. It vanishes at 
the sharp point since = 0. Therefore, with Zi = xi + zwi and Z 2 = X 2 + j(jJ 2 , |'2q/(zi, Z 2 )| < 

|zi — z/2 — 1/2| \z 2 — j 12 — 1/2|. This implies 

\F{z-i,Z2)\ = \{Q^{zl))-^ exp{'KX^)Zqf^(<dj{z2))~^\ < = 1 - ( 23 ) 

\zi ij2 1 / 2 | jZ2 jj2 1 / 2 | 

And therefore it is square-integrable in the cube. □ 


In the next theorem we are going to show that this expansion is linear independent. 


Theorem 3 

If {ca} G f’^(A), u G HU and 


exp(7rza;ia;i)u exp(—7r|x — -p) exp(7rja:2W2) 4 - exp(27rza;iu;i)c^ exp(—7r|x — bp) exp(27rja;2W2) = 0 (24) 


AeA 
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or in short notation 


ue'^(x) + ^c^e;,(x) = 0 
AeA 


then ca = 0 for all A G A. 

Proof. By Proposition 2, the series in (24) converges in We apply the Zak transform and get 

- uZqe'^{x,co) 


= exp{—TTXi)&i{oJi + ixi) 


exp(27ria:iu;i)cA exp(— 7 r|x — bp) exp(27rja:2W2) 

-AeA 


6 j(a ;2 + jxi) exp{—TTxl) 


The sum converges to a function (/(x). Therefore 

ff(x) =-=-exp(7rxi)(0i(u;i+ia;i))“ [uZqe^(x,w)J (0j(u;2 + ja; 2 ))“ exp(7rx2) (25) 

We already know that ZqeQ{x,oj) = 0 for the sharp point #, whereas Zqe'^{x,oj) ^ 0. Therefore ^(x) can 
only be an element of T^(IR^, H) if and only if u is zero. This means that g(x) — 0 and therefore ca = 0 for 
A G A. □ 

Remark 5 

The last Theorem shows that the series expansion on the relaxed lattice is unique. 


Corollary 1 

We get the coefficients via 



^g/#(x,Cu) 

Zqe'^(x,uj) 


i 2 

d xd cj 


Proof. We have 


2 f f fx 

.2,1 ’ . = Firx, w) = yP exp(27ri(xipi + wi?7i))ca exp(27rj'(x2P2 + ^2m)) 

ZoGnlX, w) 

<1 ^ J A=(p.r,)eA 

Extracting the coefficients ca by using the inverse Fourier transform gives 

c\ = [ [ exp(-27ri(a;i?i + exp(-27rj(x2?2 + ^ 2 ^ 2 )) d^xd^w. 

JqJq 2:,e()(x,w) 

Using Remark 3 we obtain the desired formula. 


□ 


6 Further research 

The requirement of our signal to be an element of lTo(K.^) is a strong restriction. The main condition 
for our proof is Zqf^ to be continuous and square integrable. A possible extension for further research 
could be Sobolev spaces or modulation spaces. Another interesting topic for further research could be the 
investigation of the convergence properties of the coefficients if more points are added to the lattice. 
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